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ABSTRACT 

We present results of new, deep Suzaku X-ray observations (160 ks) of the intracluster medium 
(ICM) in Abell 1689 out to its virial radius, combined with complementary data sets of the projected 
galaxy distribution obtained from the SDSS catalog and the projected mass distribution from our 
recent comprehensive weak and strong Icnsing analysis of Subaru/ Suprime- Gam and HST/ACS ob- 
servations. Faint X-ray emission from the ICM around the virial radius (rvir ^ 15.'6) is detected at 
4.0(7 significance, thanks to low and stable X-ray background of Suzaku. The Suzaku observations 
reveal anisotropic gas temperature and entropy distributions in cluster outskirts of r5oo ^ r ^ r^ii cor- 
related with large-scale structure of galaxies in a photometric redshift slice around the cluster. The 
high temperature (~ 5.4 keV) and entropy region in the northeastern (NE) outskirts is apparently 
connected to an overdense filamentary structure of galaxies outside the cluster. The gas temperature 
and entropy profiles in the NE direction are in good agreement, out to the virial radius, with that 
expected from a recent XMM-Newton statistical study and with an accretion shock heating model of 
the ICM, respectively. To the contrary, the other outskirt regions in contact with low density void 
environments have low gas temperatures (~ 1.7keV) and entropies, deviating from hydrostatic equi- 
librium. These anisotropic ICM features associated with large-scale structure environments suggest 
that the thermalization of the ICM occurs faster along overdense filamentary structures than along 
low-density void regions. We find that the ICM density distribution is fairly isotropic, with a threc- 
dimentional density slope of —2.29 ± 0.18 in the radial range of r25oo ^ r ^ ?'500i and with — 1.24j^Q'gg 
hr rsoo r ;S rvir, which however is significantly shallower than the Navarro, Frenk, & White universal 
matter density profile in the outskirts, p oc r~^. A joint X-ray and lensing analysis shows that the 
hydrostatic mass is lower than spherical lensing one 60 — 90%) but comparable to a triaxial halo 
mass within errors, at intermediate radii of 0. 6^2500 ^ r ^ O.Srsoo. On the other hand, the hydrostatic 
mass within 0. 4^2500 is significantly biased as low as 60%, irrespective of mass models. The thermal 
gas pressure within rsoo is, at most, ^ 50-60% of the total pressure to balance fully the gravity of 
the spherical lensing mass, and ^ 30-40% around the virial radius. Although these constitute lower 
limits when one considers the possible halo triaxiality, these small relative contributions of thermal 
pressure would require additional sources of pressure, such as bulk and/or turbulent motions. 
Subject headings: galaxies: clusters: individual (Abell 1689) — X-rays: galaxies: clusters — gravita- 
tional lensing 



1. INTRODUCTION 

Clusters of galaxies are the largest self-gravitating 
systems in the universe, and X-ray observations have 
revealed characteristics of intracluster medium (ICM) 
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which carries ^ 70% mass of known baryons, such as 
temperature, ICM mass, and metals therein. However, 
the ICM emission has been detectable only to ~ 0.6 times 
the virial radius rvir (e.g. Pratt et al. 2007) due to lim- 
ited sensitivities of detectors. This means ^ 80% of an 
entire cluster volume has been unexplored in X-rays. 

According to the hierarchical structure formation sce- 
nario based on cold dark matter (CDM) paradigm, mass 
accretion flows onto clusters, in particular, along fila- 
mentary structures, are still on-going. The ICM in the 
cluster outskirts, therefore, is expected to be sensitive to 
the structure formation and cluster evolution. However, 
since the X-ray data around the virial radius of most 
clusters were unavailable, we have not yet known the 
ICM characteristics in detail, such as density, tempera- 
ture, pressure and entropy. The ICM around the virial 
radius is a fascinating frontier for X-ray observations. 

In the era of Suzaku (Mitsuda et al. 2007), detection 
of the ICM beyond 0.6 rvir has become possible thanks 
to low and stable particle background of the XIS de- 
tectors (Koyama et al. 2007). In fact, Reiprich et al. 
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(2009) measured temperature profile of Abell 2044 out to 
^200: a radius witliin which the mean cluster mass den- 
sity is 200 times the cosmic critical density, and found 
that the profile is in good agreement with hydrodynamic 
simulations. George et al. (2009) determined radial pro- 
files of density, temperature, entropy, gas fraction, and 
mass up to 2.5 /i^ Mpc of PKS 0745-191, and found 
that the temperature profile deviates from that expected 
from the ICM in hydrostatic equilibrium with a univer- 
sal mass density profile as found by Navarro, Frenk & 
White (1996, 1997, hereafter NFW profile). Bautz et al. 
(2009) also detected X-ray emission to r2oo and found 
evidence for departure from hydrostatic equilibrium at 
radii as small as r ~ 1.3 Mpc (~ rsoo). Fujita et al. 
(2008) obtained metallicity of the ICM up to the cluster 
virial radii for the first time from a link region between 
the galaxy clusters Abell 399 and Abell 401. 

A gravitational lensing study is complementary to X- 
ray measurements, because lensing observables do not 
require any assumptions on the cluster dynamical states. 
The huge cluster mass distorts shapes of background 
galaxy images due to differential deflection of light paths, 
in a coherent pattern. As a result, multiple images, arcs 
and Einstein rings caused by the strong gravitational 
field of a cluster, appear around the cluster center, which 
is so-called strong gravitational lensing effect. It is capa- 
ble of measuring mass distributions of a cluster central 
region very well. Outside the core, weak gravitational 
lensing analysis, which utilizes the coherent distortions 
in a statistical way, is a powerful tool to constrain the 
mass distribution out to the virial radius. In particular, 
the Subaru/Suprime-Cam is the best instrument to con- 
duct the weak lensing analysis (e.g. Miyazaki et al. 2002; 
Hamana et al. 2003; Broadhurst et al. 2005b; Umetsu & 
Broadhurst 2008; Umetsu et al. 2009a, b; Hamana et al. 
2009; Okabe & Umetsu 2008; Okabe et al. 2009; Oguri et 
al. 2009), thanks to its high photon-collecting power and 
high imaging quality, combined with a wide field-of-view. 

Since the strong lensing analysis is complementary to 
weak lensing one, a joint analysis enables us to deter- 
mine accurately the cluster mass profile from the center 
to the virial radius. Pioneer joint strong and weak lensing 
studies on Abell 1689 (Broadhurst et al. 2005b; Umetsu 
& Broadhurst 2008) have shown that observed lensing 
profile is well fitted by an NFW mass profile with a high 
concentration parameter which is defined as a ratio of 
the virial radius to the scaled radius. Therefore, one of 
the best targets for the study of the ICM out to cluster 
outskirts is Abell 1689. Comparisons of X-ray observ- 
ables with the well-determined lensing mass would allow 
us to conduct a powerful diagnostic of the ICM states, 
including a stringent test for hydrostatic equilibrium, for 
the first time. In addition, two mass measurements using 
X-ray and lensing analysis are of vital importance to un- 
derstand the systematic measurement bias and construct 
well-calibrated cluster mass, which will improve the ac- 
curacy of determining cosmological parameters using a 
mass function of galaxy clusters (e.g. Vikhlinin et al. 
2009; Okabe et al. 2010; Zhang et al. 2010). 

In this paper, we describe the Suzaku observations and 
data reduction in §2, and X-ray background estimation 
in §3. In §4, we show results of spatial and spectral anal- 
yses and obtain radial profiles of temperature, electron 



density, and entropy. In §5, we discuss the hydrostatic 
equilibrium around the virial radius, and compare X-ray 
observables, such as hydrostatic mass, pressure and en- 
tropy, with lensing mass derived from a joint strong and 
weak lensing analysis. We also discuss the ICM char- 
acteristics in the outskirts in light of its cosmological 
environment, such as the large-scale structure and low- 
density voids. We summarize the conclusions in §6. 

In the present work, unless otherwise stated, we use 
iJo = 71 km s-i Mpc-i {h = Hq/IOO km s'^ Mpc'^ = 
0.71), assuming a flat universe with flu = 0.27. This 
gives physical scale 1" = 3.047 kpc at the cluster redshift 
z = 0.1832 (Struble & Rood 1999). The definition of so- 
lar abundance is taken from Anders and Grevesse (1988). 
Errors are given at the 90 % confidence level. 

2. OBSERVATION AND DATA REDUCTION 

2.1. Observations 

As shown in figure 1, we performed four-pointing 
5'M2afcu observations of Abell 1689, named Offsctl (north- 
east direction), Offset2 (southeast), Offset3 (southwest), 
and Offset4 (northwest), in the end of July 2008 with 
exposure of ~ 39 ks for each pointing. The pointings 
were coordinated so that the X-ray emission centroid of 
Abell 1689 was located at one corner of each pointing 
and the entire field-of-view (FOV) of the XIS covered 
the ICM emission up to the virial radius (15f6 ^ 2.9 
Mpc, Umetsu & Broadhurst 2008). The observation log 
is summarized in table 1. All the XIS sensors (XISO, 
XISl and XIS3) were in normal clocking mode without 
window or burst options throughout the observations. 

2.2. XIS data reduction 

In the present paper, we used the XIS data created 
through version 2.2 pipe-line data processing, which is 
available from the Suzaku ftp site^°. We re-processed 
unscreened XIS event files using Suzaku software ver- 
sion 11 in HEAsoft 6.6, together with the calibration 
database (CALDB) released on 2009 January 9. We 
applied xiscoord, xisputpixelquality, xispi, and 
xistime in this order, and performed the standard event 
screenings. Bad pixels were rejected with cleansis, 
employing the option of chipcol=SEGMENT. We 
selected events with GRADE 0,2,3,4, and 6. Good 
time intervals (GTI) were determined by criteria as 
"SAA_HXD==0 && T_SAA_HXD>436 && C0R>6 && ELV>5 
&& DYE_ELV>20 && A0CU_HK_CNT3_NML_P==1 && 
Sn_DTRATE<3 && ANG_DIST<1 . 5" . Detailed proce- 
dures of the data processing and screening are the 
same as those described in The Suzaku Data Reduction 
GmdeJ^ 

Redistribution matrix files (RMFs) of the XIS were 
produced by xisrmfgen, and auxiliary response files 
(ARFs) by xissimarf gen (Ishisaki et al. 2007). As an in- 
put to the ARE generator, we prepared an X-ray surface 
brightness profile of Abell 1689 using a double-/? model 
(sum of two (3 models, King 1962). The parameters of the 
model were determined through a least chi-square fit to a 
background-subtracted and vignetting-corrected 0.5-10 
keV surface brightness profile, created from the XMM- 
Newton MOSl image of Abell 1689. Point sources were 

ftp://ftp.darts.isas.jaxa.jp/pub/suzaku/vcr2. 2/ 
http://ftooIs.gsfc.nasa.gov/docs/suzaku/anaIysis/abc/ 
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detected and removed from the image, using ewavelet 
task of the SAS software version 8.0.0 with a detec- 
tion threshold set at la. As the background, we used 
a pubhc MOSl blank-sky dataset which is available from 
the XMM-Newton website^^. Among the blank-sky files, 
we chose data with hydrogen column densities between 
1.6 X 10^° cm-2 and 2.0 x 10^" cm"^^ match the hy- 
drogen column density of Abell 1689, 1.82 x 10^° cm~^ 
(Dickey & Lockman 1990, weighted average over 1° cone 
radius). The total exposure of the selected MOSl blank- 
sky data is 125 ks. The resulting parameters of double-/? 
model are {r^, /3)=(82.2 ± 11.1 kpc, 0.97 ± 0.17) for the 
narrower component and (183.4 ± 12.3 kpc, 0.70 ± 0.01) 
for the wider component, where Tc is the core radius. Ab- 
sorption below 2 keV, caused by a carbon-dominated con- 
tamination material on the XIS optical blocking filters, 
is included in the ARFs. In the calculation of contami- 
nation by xissimarf gen, differences of the contaminant 
thickness among the XIS sensors are taken into account, 
along with its radial dependences and secular changes 
(Koyama et al. 2007). 

Non X-ray background (NXB) of each XIS sensor 
was created using xisnxbgen, which sorts spectra of 
night-Earth observations according to the geomagnetic 
cut-off-rigidity (COR) and makes an averaged spectrum 
weighted by residence times for which Suzaku stayed in 
each COR interval (Tawa et al. 2008). Integrated period 
for the NXB creation is ±150 days of the Abell 1689 
observations. We assumed systematic errors (90% confi- 
dence limit) of 6.0% and 12.5% for the NXBs (Tawa et al. 
2008) from XIS-FI (front illuminated CCD, see Koyama 
et al. 2007, XISO and XIS3 in this paper) and XIS- 
BI (backside illuminated CCD, see Koyama et al. 2007, 
XISl in this paper), respectively, and added them in 
quadrature to the corresponding statistical errors. In the 
next section, we estimate the X-ray background (XRB) 
around Abell 1689 which is crucial to detect weak signals 
from the ICM around the virial radius. 

3. ESTIMATION OF X-RAY BACKGROUND 
3.1. Selection of blank-sky observations 

In order to estimate the X-ray background (XRB), we 
selected the nearest two blank-sky Suzaku observations 
to Abell 1689; one is Q1334-0033 (Observation ID = 
702067010, hereafter Q1334) which is 6°. 33 offset from 
Abell 1689, and the other is NGC4636.GALACTIC.1 
(Observation ID = 802039010, hereafter N4636_GAL) 
which is 8°. 68 offset from Abell 1689. The exposures 
are 12 ks and 35 ks for Q1334 and NGC4636_GAL, re- 
spectively, after the event processing and screening which 
are the same as those described in § 2.2. 

The hydrogen column densities (Dickey k. Lock- 
man 1990, weighted average over 1° cone radius) of 
Abell 1689, Q1334, and N4636.GAL, are 1.82 x 10^° 
cm-2, 2.02 X 10^° cm-2, and 1.96 x 10^° cm'^, respec- 
tively. They are in good agreement with one another 
within 11%, suggesting the same level of X-ray back- 
ground in soft X-ray band. Next we checked the count 
rates obtained from ROSAT All Sky survey (RASS) 
which are available via the NASA's HEASARC web- 

http:/ /xmm2.esac.csa.int/cxtcrnal/xmm_sw_cal/background/ 
blank_sky.shtml 



sitc^-^. Since the RASS count rates in the direction of 
AbcU 1689 is contaminated by its own ICM emission, we 
took an annular region around Abell 1689 of which in- 
ner and outer radii are 0°.5 and 1°.0, respectively. The 
RASS count rates of the two blank-sky fields were defined 
as circular regions of 1° radius. As shown in table 2, the 
count rates of Q1334 and N4636_GAL match those of the 
annular Abell 1689 field within errors in all the energy 
bands, except for count rates of 1/4 keV and 3/4 keV 
bands in the N4636_GAL field, which are significantly 
higher than the corresponding rates in the Abell 1689 
field. Therefore, the cosmic X-ray background (CXB) 
of AbcU 1689 is considered to be represented well by 
both Q1334 and N4636_GAL observations, while Galac- 
tic foreground emission (GFE) should be represented 
only by the Q1334 observation. Thus, in the following 
spectral analyses, we ignored spectra of N4636_GAL be- 
low 1.5 keV. 

3.2. Spectral analysis of the blank-sky observations 

Using XSPEC12 version 12.5.0ac, we fitted NXB- 
subtracted XIS-FI (averaged over XISO and XIS3) and 
XIS-BI spectra of Q1334 and N4636.GAL simuffancously 
with a model which represents the CXB and GFE emis- 
sions. RMFs and ARFs were created for each obser- 
vations by xisrmfgen and xissimarf gen, respectively. 
The input image for the ARFs is uniform emission over a 
circular region of 20' radius. As the CXB model, we used 
a power law with a fixed photon index of 1.41 (Kushino 
et al. 2002) and a free normalization. To represent the 
GFE, we employed a model of Henley & Shehon (2008), 
which consists of three apec (Smith et al. 2001) com- 
ponents; a non-absorbed 0.08 keV component, an ab- 
sorbed 0.11 keV one, and an absorbed 0.27 keV one, 
representing the emission from the Local Bubble (LB), 
cool halo, and hot halo, respectively. The three apec 
temperatures were all fixed. The hydrogen column den- 
sity was fixed to the Galatic values of 2.02 x 10^° cm~^ 
and 1.96 x 10^° cm'^ for Q1334 and N4636_GAL, re- 
spectively. The metal abundance and redshift were also 
fixed to one solar and zero, respectively. Relative nor- 
malizations of the two halo components were tied as cool 
halo : hot halo = 1 : 0.24 (Henley & Shehon 2008). 
Thus, the GFE model has only two free parameters; nor- 
malizations of the LB and halo. As shown in figure 2, 
this CXB+GFE model reproduces the blank-sky spectra 
successfully. The fitting result is summarized in table 3. 
The 2.0-10.0 keV surface brightness of the CXB com- 
ponent is (6.3 ± 0.2 ± 0.6) X 10"^ erg s'^ cm'^ sr"! 
(90% statistical and systematic errors), where the sys- 
tematic error, 10%, refers to Appendix 1 of Nakazawa 
et al. (2009). This is consistent with the average ASCA 
result of absolute CXB intensity (Kushino et al. 2002), 
(6.38 ± 0.07 ± 1.05) X 10"* erg s'^ cm'^ sr"! (90% sta- 
tistical and systematic errors). 

4. ANALYSIS AND RESULTS 

4.1. X-ray surface brightness profile 

In order to see how far the ICM emission of Abell 1689 
is detected, we made a surface brightness profile from 
the 0.5-10 keV XIS-FI (XISO + XIS3) mosaic im- 
age (figure 1). The center was chosen to be (a, 

^'^ http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg. pi 
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(5)=(13''ll'"29.495^ -01°20'29.92") in J2000 coordinates 
as a position of X-ray emission centroid in a 0.5-10 
keV XMM-Newton MOSl image. This position is close 
to the brightest-cluster galaxy, MaxBCG J197.87292- 
01.34110 BCG, 2".4 (- 7.3 kpc) offset from the center. 
We identified point sources from the MOSl image using 
ewavelet task of the SAS software version 8.0.0 with a 
detection threshold set at 7a. As to point sources out- 
side the MOSl image, we selected visible sources in a full 
band image of ROSAT PSPC of which flux limit corre- 
sponds to ^ 2 X 10""'^"' erg s~^ in 2.0-10.0 keV. As shown 
in figure 1 (right), circular regions (1' in radius) around 
49 point sources were masked in the XIS-FI mosaic im- 
age (their positions and flux are summarized in table 9 
in Appendix B). 

In figure 3, we show 0.5-2 keV and 2-10 keV raw 
(background inclusive, but point-source excluded) sur- 
face brightness profiles as black crosses. We created cor- 
responding 0.5-2 keV and 2-10 keV NXB profiles from 
an NXB mosaic image with xisnxbgen. The 0.5-2 keV 
and 2-10 keV NXB profiles are shown in magenta in fig- 
ure 3. Then, we obtained an XRB mosaic image. Based 
on the CXB and GFE models of AbeU 1689 described in 
§ 3.2, we simulated XRB (CXB+GFE) images of the four 
offset observations using xissim with exposures 10 times 
longer than those of actual observations. The 0.5-2 keV 
and 2-10 keV XRB profiles are shown in green in figure 3. 
The CXB level of Abell 1689 is not the same as that of 
the blank-sky fields because point sources are excluded 
in Abell 1689. Therefore, intensity of the CXB model 
when the flux is integrated to our flux limit, 2 x 10^^^ 
erg s~^ in 2.0-10.0 keV, was calculated from equation 
6 of Kushino et al. (2002), and used in the simulation. 
The resuhing CXB intensity of AbeU 1689 is 65% of the 
absolute CXB intensity. 

Finally, we obtained a background-subtracted profile 
(black — magenta — green in figure 3), as shown in red in 
figure 3. In the error bars in this background-subtracted 
profile (red crosses), Icr uncertainties of the XRB and 
NXB are added in quadrature to the corresponding sta- 
tistical 1(7 errors. We adopted 3.6% as the Icr uncertainty 
of NXB (Tawa et al. 2008), while that of XRB was cal- 
culated from la statistical errors of normalizations in 
the CXB and GFE models, determined from the spec- 
tral fitting of blank sky observations (§ 3.2). These are 
2.1% and 15.8% for the CXB and the GFE, respectively, 
corresponding to 8.5% of 0.5-10 keV XRB fiux. As a 
result, the background-subtracted profiles both in 0.5-2 
keV and 2-10 keV bands extend up to the virial radius 
(15f6 ^ 2.9 Mpc). However, as we see in § 4.2, more care- 
ful treatment of uncertainties related to a rather broad 
Suzaku point spread function (PSF), ~ 2' in half power 
diameter (hereafter HPD), is needed to evaluate the sig- 
nificance of the signal around the virial radius. 

4.2. Significance of the signal 

The signals around the virial radius found in the X-ray 
surface brightness profile (red crosses in figure 3) have 
to be tested for systematic uncertainties caused by (1) 
contaminated signals from inner ICM emission, and (2) 
residual signals from the excluded point sources. 

As for the item (1), significant amount of signals from 
inner regions of the ICM is expected to contaminate 
the signal around the virial radius, since the PSF of 



Suzaku is 2' in HPD (Serlemitsos et al. 2007). How- 
ever, as described in detail in Appendix A, the simu- 
lated double-beta profile determined from the 0.5-10 keV 
XMM-Newton MOSl image (§ 2.2) is consistent with 
the observed profile, and the contaminated signals which 
originate inside 10' contributes only 22% of signals out- 
side 10'. Therefore, the signal observed around the virial 
radius cannot be explained by the contaminated signals 
from the inner ICM emission. 

The item (2) is also due to the 2' HPD of Suzaku. 
Although we excluded circular regions around the point 
sources, half of signals from point sources is expected to 
escape from the excluded circular regions of 1' radius. As 
shown in figure 4, we simulated this residual signals from 
point sources. Details of the simulation are described 
in Appendix B. In 0.5-2 keV band, the contribution of 
the point sources in the 10f0-18.'0 region is only 26% of 
the observed profile. Meanwhile, the observed 2-10 keV 
profile can be explained well by the point sources. 

By combining the contributions of (1) and (2), we con- 
clude that in 0.5-2 keV band, these effects explain at 
most 48% of the signal around the virial radius (lO.'O- 
18.'0), and the detection is significant at A.Oa level. On 
the other hand, in 2-10 keV range, the signal is not signif- 
icant. This indicates that the detected ICM temperature 
around the virial radius might be low, compared with the 
central temperature of ^ 9 keV (Snowden et al. 2008). In 
table 4, we show signal significances when X-ray surface 
brightness profiles are evaluated separately for the four 
pointing directions. In the Offset 1 direction (northeast), 
the signal in 2-10 keV band is also detected with 2.1a 
level, suggesting harder spectrum in this direction. 

4.3. Spectral analysis 

Since we surely detected the ICM emission out to the 
virial radius from the analysis of X-ray surface bright- 
ness, we proceeded to spectral analysis. For each point- 
ing (Offsetl through Offset4), we divided the XIS image 
(figure 1) to five concentric annular regions centered on 
the X-ray emission centroid (§ 4.1); their inner and outer 
radn are 0f0-2.'0, 2f0-4.'0, 4f0-6.'0, 6f0-10.'0, and lO.'O- 
18.'0. We have four azimuthal regions at the same dis- 
tance from the center. The 49 point sources were masked 
by circles of 1' radius as we did in obtaining the X-ray 
surface brightness profile (§ 4.1). After extracting spec- 
tra of the XISO, XISl, and XIS3 from each region, we 
averaged XIS-FI (XISO and XIS3) spectra. In the OfO- 
2.'0 and 2.'0-4.'0 regions, aU spectra of Offsetl, an XIS-FI 
spectrum of Offset2, and an XIS-BI spectrum of Offset4 
are unavailable because regions irradiated by the calibra- 
tion source are removed from the analysis. 

As shown in figure 5, we used XSPEC12 version 
12.5.0ac and included the XIS-FI and XIS-BI spectra 
from the regions at the same distance, and correspond- 
ing NXB and response (RMF and ARE) files created 
through the procedure described in § 2.2. In 10.'0-18.'0 
analysis, for example, we included four XIS-FI source 
spectra and four XIS-BI source spectra. In order to de- 
termine the XRB, we also included the spectra of blank- 
sky fields (figure 2) with their NXB and response files 
(§ 3.2). Then, the NXB-subtracted source spectra were 
fitted with an absorbed thin-thermal emission model rep- 
resented by phabs x apec, added to the XRB model 
(§ 3.2), while the blank-sky spectra were fitted with the 
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XRB model simultaneously with the source spectra. Rel- 
ative CXB normalization of Abell 1689 to the blank- 
fields was fixed to 65% (§ 4.1). Free parameters in the 
apec model are temperature, metal abundance, and nor- 
malization. The hydrogen column density was fixed to 
the Galactic value of 1.82 x lO^*' cm^^ and redshift was 
also fixed to 0.1832. These free parameters for different 
pointings (Offsetl through Offset4) were determined sep- 
arately when we saw azimuthal differences, or tied when 
we obtained azimuthal averages. Relative normalization 
between the XIS-FI and XIS-BI was left free and its de- 
viations from unity were within ±9% in all the regions. 
The fitting result is summarized in table 5. Although the 
X-ray background model is determined separately for the 
individual annular regions, the normalizations of CXB, 
LB, and Halo are consistent among the regions within 
errors. In this spectral analysis, we did not explicitly 
include the two uncertainties described in § 4.2, namely, 
contaminated signals from inner ICM emission and resid- 
ual signals from the excluded point sources. We estimate 
their effect on the temperature and electron number den- 
sity of Abell 1689 in the following section, obtaining ra- 
dial profiles of these ICM parameters. 

4.4. Radial profiles 
TEMPERATURE PROFILE 

The top panel of figure 6 shows radial profiles of gas 
temperature, obtained from the spectral analysis (§ 4.3). 
In figure 6 (and following figures of any radial profiles), 
the center of each radial bin is defined as emission- 
weighted center derived from the X-ray surface bright- 
ness obtained in § 4.1; l.'l, 2'.7, 4'.8, 7.'8, and 13.'4 for the 
five annular regions. Although the temperature profile 
is a projected one obtained by the two-dimensional spec- 
tral analysis, it coincides well with a three-dimensional 
temperature profile obtained with Chandra data (Peng 
et al. 2009) in the overlapping region. The temperature 
is ~ 10 keV within 4.'0 730 kpc), and it gradually de- 
creases toward the outskirts down to ^ 2 keV around the 
virial radius (rvir ^ 15f6 ^ 2.9 Mpc). The decrease of 
temperature is independent of the azimuthal directions 
up to ~ rsoo ~ 8.' 7, at which the mean interior density is 
500 times the critical mass density of the universe. Note 
that the overdensity radius rsoo is determined by lensing 
analysis, as is the case for rvir (see §5.2 and eq. (7)). 
The logarithmic slopes of temperature between the third 
and forth bins, corresponding to r25oo ^ r ^ t^qq are 7 = 
-dlogT/dlogr = 0.87 ±0.25, 0.71 ±0.30, & 0.87±0.15 
for all, Offsetl and Offset234 directions, respectively, 
which are consistent within errors with each other. 

The temperature variations in direction becomes sig- 
nificant in the outskirts of r5oo ^ r ^ rvir ■ the tempera- 
tures for Offsetl (~ 5.4 keV) and Offset234 (~ 1.7 keV) 
are about half and one-sixth of a mean temperature 10.4 
keV, which is an emission- weighted temperature of 2'.0~ 
4f0 and 4f0-6f0 regions, respectively. This temperature 
variations can be seen directly in spectra. In figure 7, 
we show the spectra of Offsetl and Offset2 in 10f0-18.'0 
region. The signal of Offsetl is clearly detected from 
both XIS-FI and XIS-BI in 2-5 keV range, whereas the 
signal of Offset2 is comparable to the background in the 
same energy range. This means that the XIS spectra 
of Offsetl is harder than those of Offset2, resulting in 



the higher temperature in the Offsetl direction. The 
logarithmic slopes of temperature in the outskirts cor- 
responding to r5oo ^ r a fvir are 7 = — dlogT/dlogr — 
1.45 ± 0.46, 0.18 ± 0.41 & 1.92 ± 0.24 for aU, Offsetl and 
Offset234 directions, respectively. The slope in Offset234 
and all directions are surprisingly steep. 

Contaminated signals from the inner ICM emission due 
to the Suzaku PSF may affect the temperature determi- 
nation. However, in PKS 0745-191 case (George et al. 
2009), the effect is at most AT = -0.8 keV even in a 
central region where > 30% of photons come from the 
other regions. In the other regions where photon mixing 
is < 20%, the temperature uncertainties are < 0.1 keV. 
Therefore, in our case of Abell 1689, we expect that the 
PSF effect is at most ~ 0.1 keV in the cluster outskirts 
outside 6', which is much smaller than the statistical er- 
rors. The effect inside 6' would be also negligible because 
our temperature profile is consistent with the Chandra 
result (Peng et al. 2009) as described above. In order 
to see the other effect, the residual signals from the ex- 
cluded point sources in 10f0-18.'0, we added a power-law 
of photon index 1.4 into the fit with its normalization 
fixed to the summed escaped fiux of the 49 excluded point 
sources. This corresponds to increasing the CXB level by 
^ 15%. As a result, the effect was as small as AT ^ —0.3 
keV in all the azimuthal directions, which is also smaller 
than the statistical errors. 

We compared our temperature profiles with that based 
on the scaled temperatures derived from a sample of 15 
clusters by Pratt et al. (2007). We used the mean tem- 
perature 10.4 keV to calculate the scaled profile. The 
temperature profile in the Offsetl direction is in good 
agreement with the scaled temperature one (dashed line 
in figure 6(a)), while the temperatures of the other re- 
gions (Offset2 - Offset4) around the virial radius are lower 
than the scaled value by a factor of ~ 2. 

DENSITY PROFILE 

The electron number density profile was calculated 
from the normalization parameter of apec model, defined 
as 

10"" f 

"^"^ = M^KTTzjr J 

where Da is the angular size distance to the source in 
units of cm, rie is the electron density in units of cm~^, 
and njj is the hydrogen density in units of cm"'^. Since 
metal abundance was poorly constrained from our fitting 
result (table 5), we used nn/uc = 0.852 assuming metal 
abundance to be 0.3 solar. We calculated a geometrical 
volume which contributes to each two-dimensional region 
assuming that Abell 1689 is spherically symmetric, and 
obtained projected volume-averaged density for the re- 
gion. Then, the projected density profile was deprojected 
by the following equation 

y "c.dcpro^^i = "-c.proj^ ' (2) 

where 

Tic.dcpro is the deprojected electron density, ?T'c,proj 
is the projected electron density, and Vi is the volume 
contributing to i-th annular region. We assumed that no 
emission exists outside the 10'0-18.'0 region and hence 
'T-o.dopro = ^o.proj in the 10f0-18f0 rcgiou. The result is 
summarized in table 6. 



6 



Kawaharada et al. 



The electron number density profile (middle panel of 
figure 6) decreases down to 0(10"'^) cm"'^ around the 
virial radius. The density profiles arc in good agree- 
ment with a generalized /3 model (Vikhlinin et al. 1999) 
obtained by fitting Chandra data within 9'.0 (Peng et 
al. 2009, Model 3 in table 9). In contrast to the tem- 
perature profile, the electron densities of the four off- 
set directions are consistent with one another within er- 
rors in any circular annulus. The slope of the number 
density is —2.29 ± 0.18 in the range of r25oo ^ r ^ r^oo- 
Our result coincides with previous studies around rsoo 
: oc r-^-^ for the REFLEX-DXL clusters (Zhang et 
al. 2006, 2007), Ue oc r-^-^io.i foj. ^j^g LoCuSS sample 
(Zhang et al. 2008), oc r-^-80±o.46 ^-j^^, REXCESS 
sample (Croston et al. 2008), and rig oc r~^-^^ for 13 
clusters observed with Chandra (Vikhlinin et al. 2006). 
A consistent result, Ue oc 7.-2-27±o.i2^ j^g^g i^ggjj g^jgQ qI^. 
tained with Suzaku in Abell 1795 around rsoo (Bautz et 
al. 2009). In the outer annulus of rsoo ^ r rvir, the 
density slope becomes flatter to oc 7.-i-24±o.23^ r^j^^ 
same tendency of flattening is seen in a electron density 
profile of PKS 0745-191 around r > 2 Mpc (George et al. 
2009). Excess signals in a X-ray surface brightness pro- 
file above a (3 model outside 2 Mpc seen in Abell 1795 
(Bautz et al. 2009) also indicates the fiattening of slope in 
the cluster outskirts. The electron density profile outside 
T ^ fboo is significantly shallower than asymptotic mat- 
ter density slope p oc r~^ of an NEW profile (Navarro, 
Frenk, & White 1996, 1997), defined by 



^ (r/r,)(l + r/r,)2' ' 

where ps is the central density parameter and r,, is the 
scale radius. 

Since the electron density profile is consistent with the 
Chandra result within 9f0 (Peng et al. 2009), we can say 
that the systematic effects of the Suzaku PSF are safely 
within the statistical errors in this region. In our lO.'O- 
18.'0, however, 48% of the signal can leak from the other 
regions or point sources as discussed in § 4.2. This frac- 
tion corresponds to 28% in terms of the electon density, 
which is comparable to the statistical error in this region 
(see table 6). The steepest slope including this system- 
atic effect is —1.8, and is still much shallower than the 
NEW profile. 

Cavaliere et al. (2009) analytically constructed a phys- 
ical model of the ICM in hydrostatic equilibrium with 
the dark matter potential, in the context of dark-matter 
Jeans equilibria (Lapi & Cavaliere 2009) constrained 
by the dark matter entropy in a single power-law form 
K{r) oc r", as motivated by N-body simulations of cold 
dark matter halos (Taylor & Navarro 2001; Faltenbacher 
et al. 2007). Their equilibrium model predicts a shallow 
gas density slope of (dlnn/dlnr ^ —2.2) at the entropy 
slope of ~ 1.1 in a slow and constant accretion phase, 
in good agreement with our observed gas density slope 
of —2.29 ± 0.18 in the range of r25oo ^ ''soo- Beyond 
r5oo of Abell 1689, we calculated the ratio of the dark 
matter circular velocity to the gas sound velocity from 
the observed gas density and entropy slopes, by following 
their model (equation 7 in Cavaliere et al. 2009). The 
resulting ratio is ~ 1 — 2, which conflicts with our ob- 
servations that the dark matter velocity is much higher 



than the sound velocity. As we discuss in § 5, the ICM 
in most regions at the outermost annulus is likely not 
to be in hydrostatic equilibrium. Therefore, some mod- 
ifications to the model would be needed to describe our 
results outside rsoo- 

ENTROPY PROFILE 

The bottom panel in figure 6 shows the entropy profile 
calculated as 



using the temperature and electron density profiles ob- 
tained above. 

The entropies in all directions increase from the center 
to the radius of lO.'O. The entropy distribution in the last 
radial bin is anisotropic, similarly to the temperature 
profiles. The radial dependence of the Offset 1 entropy 
obeys roughly r^-^ (dashed line in figure 6(c)), which is 
predicted by models of accretion shock heating (Tozzi & 
Norman 2001; Ponman et al. 2003), while the entropies 
of the other regions (Offset2-Offset4) in the last radial 
bin are lower than those in their 6.'0-10f0 annuli. 

5. DISCUSSION 

5.1. Thermodynamics in the Outskirts 

In this section we discuss the distributions of gas tem- 
perature and entropy shown in figure 6, particularly fo- 
cusing on their anisotropic distributions in the last bin 

(^500 ^ r < fvir). 

In the framework of the hierarchical clustering scenario 
based on CDM paradigm, mass aggregation processes 
onto clusters, in particular along the large-scale filamen- 
tary structures, are still on-going. Continuous mass ac- 
cretion flows along filamentary structures, within which 
clusters are embedded, sometimes trigger internal shocks 
in gas within clusters. It converts most of the tremen- 
dous amounts of kinetic energy into the thermal energy 
to heat the ICM. Therefore, the accretion fiows are con- 
sidered to play an important role in cluster evolution as 
well as in the ICM thermodynamics. If the shock heating 
dominates in the cluster outskirts, the entropy generated 
in the shock process becomes higher than those in the 
pre-shock regions. On the other hand, if the gas is adi- 
abatically compressed during the infall with a sub-sonic 
velocity, the gas temperature increases but the entropy 
stays constant. 

According to analytical models and simulations consid- 
ering both the continuous accretion shock and adiabatic 
compression processes (e.g. Tozzi & Norman 2001; Voit 
et al. 2002; Borgani et al. 2005), the entropy profile is pre- 
dicted to increase with K oc r^-^ by the accretion shock, 
as long as the initial entropy is low, assuming that the 
ICM is in hydrostatic equilibrium and the kinetic energy 
of the infalling gas is instantly thermalized via shocks. 
This radial dependence has been generally confirmed in 
observed profiles (e.g. Ponman et al. 2003; Pratt et al. 
2006) within rsoo. Interestingly, our entropy profile in 
the Offsetl direction coincides very well with the model 
prediction up to the virial radius. It indicates a possibil- 
ity that the ICM in the Offsetl rsoo ^ r ^ rvir region is 
thermalized and close to the hydrostatic equilibrium. In- 
deed, the outskirts temperature in the Offsetl direction 
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agrees with the predictions of the analytical models (e.g. 
Komatsu & Seljak 2001; Tozzi & Norman 2001; Ostriker 
et ah 2005) in which, under the hydrostatic equiUbrium 
assumption, the temperature around the virial radius de- 
chnes at most 50 percent of the inner values. 

Contrary to the Offset 1, the sharp decline in both tem- 
perature and entropy in the Offset2-Offset4 regions indi- 
cates that the thermal pressure is not sufficient to bal- 
ance the total gravity of the cluster. If the gas is not 
completely thermalized, the gas temperature should be 
lower than the value needed to maintain hydrostatic equi- 
librium. In this case, the accreted gas retains some frac- 
tion of the total energy in bulk motion, and the bulk 
pressure partly supports the gravity. In fact, recent nu- 
merical simulations (e.g. Vazza et al. 2009) have shown 
that the kinetic energy of bulk motion carries ^ 30% of 
the total energy around the virial radius. In addition to 
the bulk pressure, pressure of turbulent motions in the 
ICM may also contribute to balance the total gravity 
(e.g. Nagai, Vikhlinin, & Kravtsov 2007; Piffaretti & 
Valdarnini 2008; Jeltema et al. 2008; Lau, Kravtsov, & 
Nagai 2009; Fang et al. 2009; Vazza et al. 2009). 

In order to study the validity for the assumption of hy- 
drostatic equilibrium, we measured the hydrostatic equi- 
librium masses (H.E. mass) for the all, Offsetl and Off- 
set234 regions. We here assume that the mass distribu- 
tion is spherically symmetric and the ICM is in hydro- 
static equilibrium at any radii. The hydrostatic mass, 
-^^H.E., is calculated from the parametric density and 
temperature profiles obtained by fitting data points in 
each azimuthal region (all, Offsetl, & Offset234) with 
models (Appendix C and table 7), as follows 

Pg dr 

Pg = —UekBT, 

A* 

where Pg is the gas pressure and pg = fiempUe is the 
gas mass density. The errors (68% CL uncertainty) are 
calculated with Monte Carlo simulations because model 
parameters are correlated with each other. The resulting 
mass is shown in figure 8. The cumulative mass in the 
Offsetl direction increases with radius, while the ones in 
all and Offset234 regions unphysically decrease outside 
~ 7f0~8.'0. This unphysical decrease in the cumulative 
mass is due to the sharp temperature drops in the Off- 
set234 region. Here, in the parametric density model 
(equation CI), we adopted a modified (3 model to re- 
produce the density flatness in outskirts (r > lO.'O). The 
ratio of core radii rc.2/'i'c.i in equation CI was fixed to be 
3, since the fiattening factor in which rc,2 is included was 
not constraind well owing to few data points. However, 
even by choosing the ratio A{rc^2/fc,i) = il, the radii, 
at which the masses have maximum values, are changed 
only by +7%, -3% for Offset234 and +12%, -6% for ah 
regions, respectivly. This is because the curvature of the 
parametric density profile is insensitive to a choice of the 
ratio of core radii. Thus, the hydrostatic equilibrium we 
assumed is inadequate to describe the ICM in the out- 
skirts except the Offsetl direction, that is, most of the 
ICM in the cluster outskirts is far from hydrostatic equi- 
librium. 

One alternative proposal to the bulk motion or turbu- 



lent motion as a cause of the non- hydrostatic equilibrium 
would be convective instability in the cluster outskirts. 
We investigated the possibility of convective instability 
in the radial direction. The convective motion in Off- 
sct234 is unstable outside ^ 9f0 but the time scale of the 
growing mode is comparable to the age of the universe 
at cluster redshift z = 0.1832. Therefore, the convective 
instability is negligible for the unphysical decrease in the 
cumulative mass. Another possible proposal is that ions 
have a higher temperature than that of electrons (Fox 
& Locb 1997; Ettori & Fabian 1998; Takizawa 1998), 
and that the thermal pressure of ions supports the clus- 
ter mass. We computed the thermal equilibration time 
between electrons and ions through the Coulomb inter- 
action. In the cluster outskirts of Abell 1689, the time 
scale is given by 

(e.g. Spitzer 1962; Takizawa 1999; Akahori & Yoshikawa 
2009). The Coulomb interaction would provide us with 
the maximum time scale of thermal equilibrium, because 
there might be other processes, such as plasma insta- 
bility, to facilitate the interaction between electrons and 
ions. Since the resulting time is much shorter than the 
dynamical time scale, it is difficult to explain the un- 
physical decrease. Therefore, the bulk and/or turbulent 
motions would be the main source (s) of additional pres- 
sure needed in the Offset234 directions. In §5.2, we shall 
discuss on this point again based on a joint X-ray and 
lensing analysis. 

What makes the anisotropic distributions of gas tem- 
perature and entropy in the outskirts? The internal 
shocks associated with mass accretion flows would heat 
the ICM. Numerical simulations (e.g. Ryu et al. 2003; 
Molnar et al. 2009) have shown that, the low-density 
gas in low-density regions, so-called voids, accretes onto 
the vicinity of clusters by the gravitational pull, with 
an order of thousands of kilometers per second. As a 
result, the internal shocks form outside virial radius at 
which gas density sharply changes. It is also referred to 
as virial shocks (Ryu et al. 2003). Meanwhile, the ac- 
cretion flow along filamentary structures forms internal 
shocks inside the virial radius. These shocks migrate rel- 
atively deep into a cluster potential well, because a large 
amount of matter, such as gas, galaxies and groups, ac- 
cretes from filaments (Ryu et al. 2003). Therefore, the 
thermalization in the outskirts along the direction of fil- 
amentary structure takes place faster. Based on this ac- 
cretion shock scenario, the observed anisotropic distri- 
butions of gas temperature and entropy in the outskirts 
inside virial radius would be associated with large-scale 
structures. Another possibility for the anisotropy is that 
infalling galaxies, especially massive galaxies, might be 
reservoirs of cold gas. The cold and dense gas, which is 
bound by the dark matter potential of galaxies, would 
decrease the emission-weighted temperature in the out- 
skirts. Although we excluded point sources in the spec- 
tral study, we cannot rule out this possibility only from 
the X-ray results. If this is the case, member galaxies 
around the virial radius (10.'0-18f0) in the Offsct234 di- 
rections would be more numerous than those in the Off- 
setl direction. 
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In order to search the large-scale structures or the 
anisotropic galaxy distribution around the virial radius, 
we made maps of galaxy number density using the SDSS 
DR7 catalogue (Abazajian et al. 2009) from SDSS 
CasJobs site-'^"'. We first looked into bright red-sequence 
galaxies which are expected to host cold gas. We selected 
red-sequence galaxies by criteria of \ {g' — i') — (— 0.048i' + 
2.556)1 < 0.152 and i' < 21 ABmag, where we used psf- 
Mag for magnitude and modelcolor for color. We made 
projected galaxy distribution, smoothed with a Gaussian 
kernel, w = exp(-rVr2), of FWHM= 2{ln2)^^^rg = 
l'.67. The resulting map is shown in the left panel of 
figure 9. We also looked into the same map derived from 
Subaru data. It is consistent with the SDSS map, al- 
though a part of the outer region is not available in the 
Subaru data due to the limitation of its field-of-view. 
In figure 9 (left), the central part of density map is 
clearly elongated in the north-south direction. Umetsu & 
Broadhurst (2008) have found the same elongation in a 
two-dimensional mass distribution. However, we could 
find neither an apparent feature nor a significant dif- 
ference between the projected distributions in the Off- 
set 1 and the other directions in the cluster outskirts of 
lO'.O < r < 18'.0. Therefore, the temperature and en- 
tropy anisotropics in the outskirts are not due to the 
cold gas in galaxies. 

Next, we made a larger map to investigate the pres- 
ence of large-scale structures associated with Abell 1689. 
Lemze et al. (2009) has shown that the maximum line- 
of-sight velocity of infall galaxies is ^ 4000 km/s. There- 
fore, taking into account the uncertainty of rcdshift due 
to line-of-sight velocities, we selected galaxies in the 
range of \z — Zc\ < Sz = cri,,max(l + ^c)/c — 0.0158, 
where z^ is the cluster redshift, z is a photometric red- 
shift, av,max = 4000 km/s, and c is the light veloc- 
ity. The mean photometric redshift for our sample is 
z = 0.18320 ± 0.00004 ± 0.00035, where the first error 
is the statistical error and the second one is the system- 
atic error due to photometric uncertainty of each galaxy. 
Our sample of galaxies thus statistically represents a slice 
around Abell 1689 in redshift space. 

The resulting map of galaxies smoothed with a Gaus- 
sian kernel of FWIIM=:20f0 is shown in the right panel 
of figure 9. In the Offsetl direction, a filamentary over- 
density region outside the virial radius is found, where 
the galaxy number density is more than 1.5 times higher 
than those in the other directions. In the Offset4 direc- 
tion outside the FOV of XIS, a narrow sheet of galaxies is 
also found, which is connected to a northwest overdensity 
region. We note that the projected elongated direction 
of cluster mass and galaxies (left panel of figure 9) does 
not coincide with the filamentary direction (right panel), 
suggesting that a mass structure seen in the central re- 
gion of a cluster is not necessarily connected directly to 
a filament. Our result does not change even when we 
choose twice or half Sz in the galaxy selection. We also 
tried to investigate the line-of-sight filamentary structure 
by a Monte -Carlo simulation computing the redshift dis- 
tribution of galaxies with photometric errors. However, 
we could not obtain reliable results within a rcdshift res- 
olution smaller than ^ 10 Mpc. We could not see the 

http://casjobs.sdss.org/ 



line-of-sight structure from the spectral data of SDSS ei- 
ther, because available number of galaxies is too small. 

The high temperature and entropy region in the out- 
skirts is clearly correlated with the galaxy overdensity 
region associated with the large-scale structure outside 
Abell 1689, as shown in figure 10. This indicates that 
the ICM in the outskirts is significantly affected by sur- 
rounding environments of galaxy clusters, such as the 
filamentary structures and the low density void regions. 
The large-scale structure would play an important role in 
the thermalization process of the ICM in the outskirts. In 
particular, our result suggests that the thermalization in 
the outskirts along the filamentary structure takes place 
faster than that in the void region. 

5.2. A Joint X-ray and Lensing Analysis 

In this subsection, we carry out a joint X-ray and lens- 
ing analysis, incorporating Subaru/ Suprime- Cam and 
HST/ACS data. Abeh 1689 has been the focus of inten- 
sive lensing studies in recent years (e.g, Broadhurst et 
al. 2005a, b; Limousin et al. 2007; Umetsu & Broadhurst 
2008; Corless et al. 2009). It has been shown by Broad- 
hurst et al. (2005b) and Umetsu & Broadhurst (2008) 
that joint lensing profiles of Abell 1689 obtained from 
their ACS and Subaru data with sufficient quality are 
consistent with a continuously steepening density profile 
over a wide range of radii, r = 10 — 2000kpcft.~^, well 
described by the general NEW profile (Navarro, Frenk, 
& White 1996; Navarro, Frenk, & White 1997), whereas 
the singular isothermal sphere model is strongly disfa- 
vored. They also revealed that the concentration param- 
eter of the NFW profile, namely the ratio of the virial 
radius to the scale radius, Cvir = fvii/^'sj is much higher 
than predicted by cosmological iV-body simulations (e.g.. 
Bullock et al. 2001; Neto et al. 2007). We note that the 
virial radius, Tvir, within which the mean interior density 
is Avir — 110 (Nakamura & Suto 1997) times the critical 
mass density, pcr{z), at a cluster redshift is given by 

Mvir = ^7rAvirPcr(2)?'vir- (7) 

In the present paper, the full lensing constraints derived 
from the joint ACS and Subaru data allow us to com- 
pare, for the first time. X-ray observations with the total 
mass profile for the entire cluster, from the cluster cen- 
ter to the virial radius. The lensing analysis is free from 
any assumptions about the dynamical state of the clus- 
ter, so that a joint X-ray and lensing analysis provides a 
powerful diagnostic of the ICM state and any systematic 
offsets between the two mass determination methods. 

Here we first summarize our lensing work on Abell 1689 
before presenting the results from our joint analy- 
sis. Umetsu & Broadhurst (2008) combined HST/ACS 
strong lensing data with Subaru weak lensing distor- 
tion and magnification data in a two-dimensional anal- 
ysis to reconstruct the projected mass profile. Their 
full lensing method, assuming the spherical symme- 
try, yields best-fit NFW model parameters of A/vir = 
l-47lo:33 X IQi^Mo/i-i and c^r = 12.7 ± 2.9 (includ- 
ing both statistical and systematic uncertainties; see 
also Lemze et al. (2009)), which properly reproduce 
the observed Einstein radius of 9e = 45" for Zg ~ 1 
(Broadhurst et al. 2005a). Here we deproject the two- 
dimensional mass profile of Umetsu & Broadhurst (2008) 
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and obtain a non-parametric M3D profile simply assum- 
ing spherical symmetry (Broadhurst & Barkana 2008; 
Umetsu et al. 2009b). This method is based on the 
fact that the surface- mass density Em(-R) is related to 
the three-dimensional mass density p(r) by an Abel inte- 
gral transform; or cquivalently, one finds that the three- 
dimensional mass M3d(< r) out to spherical radius r is 
written in terms of (i?) as 

M3d(< ^) = dRm:^{R)~A dRRf ^m{R), 

(8) 

where /(x) = (x^ - l)-i/2 - tmi-^{x^ - (Broad- 
hurst & Barkana 2008). Wc propagate errors on S„i(i?) 
using Monte-Carlo techniques taking into account the 
error covariance matrix of the full lensing constraints. 
This deprojection method allows us to derive in a 
non-parametric way three-dimensional virial quantities 
(rvirjMvir) and values of Ma = M3d(< ^a) within a 
sphere of a fixed mean interior overdensity A with re- 
spect to the critical density of the universe at the cluster 
redshift. The resulting parameters of the non-parametric 
spherical, as well as NFW, models are summarized in ta- 
ble 8. 

MASS COMPARISON 

Here we compare our X-ray hydrostatic mass measure- 
ments (in all regions; §5.1) with the spherical-lensing 
and NFW mass profiles in figure 11. As discussed in 
§5.1, the hydrostatic equilibrium assumption is invalid 
in the cluster outskirts and hence, a mass comparison 
is limited within the radius ^ 7', inside which the hy- 
drostatic mass increases with radius. The hydrostatic 
mass is systematically lower than the lensing mass at all 
radii. In particular, the hydrostatic mass in the central 
region (r ^ 2f0) is significantly lower than the lensing 
one (AfH.E.Z-^^icns ~ 60%). Since our hydrostatic mass 
estimates agree with the latest Chandra results (Peng et 
al. 2009) in the overlapping region, this would not be 
due to the low angular resolution of the Suzaku satel- 
lite. At intermediate radii of 3f0 r ;S 7.'0, the hy- 
drostatic to lensing mass ratio, Mh.e. /Miens, is in a 
range of ~ 60 — 90%. Some numerical simulations (e.g. 
Nagai, Vikhhnin, & Kravtsov 2007; Piffaretti & Val- 
darnini 2008; Fang et al. 2009) have shown that X-ray 
hydrostatic-equilibrium mass estimates would underes- 
timate the actual cluster mass, because the turbulent 
pressure and the bulk/rotational kinetic pressure par- 
tially support the gravity of the cluster. For instance, 
Nagai, Vikhlinin, & Kravtsov (2007) found that the hy- 
drostatic mass is biased low, and is underestimated on 
average by 0.121±0.136 and 0.163±0.095 at overdensities 
A = 2500 and 500, respectively, for a sample of 12 sim- 
ulated clusters at redshift z = 0. Piffaretti & Valdarnini 
(2008) found that the hydrostatic mass derived from the 
/3-model gas density profile and spectroscopic-like tem- 
perature model, similarly to our modeling, is biased low 
by 12 ± 11% and 27 ± 11% at overdensities of A = 2500 
and 500, respectively. They found that, on average, the 
systematic underestimate bias of Mh.e. (with respect to 
the actual cluster mass) increases with decreasing A, or 
increasing the pivot radius ta. 

We find that the degree of systematic bias in the X- 



ray mass estimate, Mjj.e. /Miens, is in rough agreement 
with the results from numerical simulations, except for 
the central region. However, such a comparison with 
simulated clusters in the central region would be practi- 
cally difficult owing to over-cooling problems in cosmo- 
logical hydrodynamical simulations, failing to reproduce 
the observed ICM features. If the ICM is not in ther- 
mal balance with the gravity and an additional pressure 
support is provided by the gas bulk/turbulent motions, 
as suggested by simulations, the X-ray micro-calorimeter 
instrument (SXS) on the next Japanese X-ray satellite 
ASTRO-H (Takahashi et al. 2008) will directly detect 
them and measure their contributions to the pressure 
balance. 

Another possible cause of the mass discrepancy is the 
triaxial structure of clusters. If the major axis of a tri- 
axial halo is pointed to the observer, then this orien- 
tation boosts the projected surface mass density S],„(0) 
and hence the lensing signal. As a result, halo triax- 
iality and orientation bias can affect the lensing-based 
mass and concentration measurements (e.g, Oguri et 
al. 2005). Oguri et al. (2009) demonstrated in the con- 
text of the ACDM model that clusters with large Ein- 
stein radii (^e > 20"), or the so-called superlenses, are 
likely to be highly biased, with major axes preferentially 
aligned with the line-of-sight. The level of correction for 
the orientation bias in lensing mass estimates (obtained 
a priori assuming the spherical NFW model) is derived 
from a semi-analytical representation of simulated CDM 
triaxial halos by Oguri et al. (2005). Oguri et al. (2005) 
fitted the two-dimensional mass map of Umetsu & Broad- 
hurst (2008) with a triaxial halo model and obtained 
Mvh. = 1.26t!]:^^ X IO^^/i-^Mq and Cvir = 16.9t?2^g (ta- 
ble 8) where the quoted errors here are statistical only 
(68.3% CL). This triaxial description of the full lensing 
constraints presented in Umetsu & Broadhurst (2008) is 
in better agreement with our X-ray results. On the other 
hand, more recently, Peng et al. (2009) demonstrated 
that a simple triaxial modeling of Chandra X-ray data 
can relax the mass discrepancy with respect to strong 
and weak lensing. They found that the total spherical 
mass and the spherically-averaged mass density are es- 
sentially unchanged under different triaxiality assump- 
tions, as previously reported by Piffaretti et al. (2003) 
and Gavazzi (2005). By assuming a prolate symmetry 
of the gas structure with axis ratio 0.7, they found that 
their X-ray mass estimate agrees with those of strong 
and weak lensing within 1% (— lo") and 25% (-flcr), re- 
spectively. 

Using the triaxial halo model of Oguri et al. (2005), 
the X-ray to lensing mass ratio becomes consistent with 
unity at larger radii within the errors; however, it is still 
much lower than unity in the central region due to the 
inferred high concentration parameter. Based on inde- 
pendent weak lensing data, Corless et al. (2009) also ob- 
tained high concentrations from their triaxial halo mod- 
eling in conjunction with the strong lensing prior on the 
Einstein radius. Due to the tight strong lensing con- 
straints that prefer a highly-concentrated mass distribu- 
tion, it is difficult to explain the mass discrepancy in the 
central region. 

GAS FRACTION 
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From our combined X-ray and Icnsing data set, we 
measure the cumulative gas mass fraction, 



where Mgas(<) and Miens (< r) are the gas and lens- 
ing masses contained within a sphere of radius r. The 
gas mass is calculated by fitting a modified (3 model 
to the electron number density profile (Appendix C). 
The errors are obtained by 1000 Monte-Carlo simula- 
tions taking into account the error covariance matrix. 
The resulting gas fraction, shown in the right panel of 
figure 11, increases toward the outskirts, but does not 
reach the cosmic mean baryon fraction (WMAP5; Ko- 
matsu et al. 2009) within the virial radius. In figure 11 
(right), we also compare our result with those from 
previous statistical studies (Zhang et al. 2010; Umetsu 
et al. 2009a): the mean gas fractions for 12 clusters 
from the Local Substructure Survey (LoCuSS) based on 
a joint XMM-Newton X-ray and Subaru/Suprime-Cam 
weak lensing analysis (Zhang et al. 2010) and those for 
4 high-mass clusters (Mvir ^ IO^^Mq h~^) from a joint 
Sunyaev- Z'eldovich effect (SZE) AMiBA and Subaru 
lensing analysis (Umetsu ct al. 2009a). We note that 
these gas fraction measurements, based on the lensing- 
derived total mass, are independent of the cluster dy- 
namical state. The mean gas fractions from the two sta- 
tistical studies are in close agreement with each other. 
On the other hand, our gas fraction measurements of 
Abell 1689 based on the joint X-ray/lensing analysis, as- 
suming spherical symmetry, are about half of the mean 
cluster values derived from Zhang et al. (2010) and 
Umetsu et al. (2009a). For the case of the triaxial mass 
model (Oguri et al. 2005), the upper bounds of the gas 
fraction come closer to, and are marginally consistent 
with, the mean values derived by statistical work. Re- 
cently Bode et al. (2009) proposed a method of modeling 
the ICM distributions in gravitational potential wells of 
clusters, including star formation, energy input, and non- 
thermal pressure support calibrated by numerical simu- 
lations and X-ray observations. In their standard model, 
the gas initially contained inside the virial radius ex- 
pands outward. Consequently, the cluster baryon frac- 
tion reaches the cosmic mean at r w 1.2rvir, which is in 
rough agreement with our result. 

PRESSURE AND ENTROPY PROFILES DERIVED 
FROM GRAVITATIONAL LENSING 

We compare gas thermal pressure obtained from the 
two X-ray observablcs, temperature and electron num- 
ber density, with a pressure calculated from the lensing 
mass profile assuming the hydrostatic equilibrium (here- 
after lensing pressure). The lensing pressure is a total 
pressure to balance fully the gravity of the cluster lens- 
ing mass (see also Molnar et al. 2010). We calculated 
the lensing pressure profile, Pions(?'), by integrating eq. 
(5) with a boundary condition that the external pres- 
sure is zero at r = 40f0 (« 20vir)- The inner pressure 
is little affected by choices of the outer boundary condi- 
tion as long as the pressure outside the virial radius is 
significantly lower than the values at inner radii. We ex- 
trapolated the observed lensing mass and density profile 
outside the field-of-view. The errors were computed by 



1000 Monte Carlo simulations. 

The lensing pressure is systematically higher than the 
gas thermal pressure at any radii (the left panel of figure 
12). Within the radius of r25oo, the gas thermal pressure 
is at most 50 - 60% of the lensing pressure. In the out- 
skirts around the virial radius, the gas pressure in Offset 1 
is comparable to lensing one, while ones in the Offset234 
and all regions is ~ 30 ~ 40%. This is inconsistent with 
a result from SZE analysis of clusters; a stacked SZE sig- 
nal profile for galaxy clusters using WMAP three-year 
data (Afshordi ct al. 2007; Atrio-Barandcla ct al. 2008) 
has shown that the gas closely follows the dark matter 
distribution out to the virial radius. 

If we assume that the mass discrepancy is fully ex- 
plained by the gradient of additional pressure supports 
such as turbulent and/or bulk pressures, they contribute 
up to ^ 40 - 50% of the lensing pressures within r5oo and 
~ 60 - 70% around the virial radius, respectively. Hydro- 
dynamic numerical simulation (Lau, Kravtsov, & Nagai 
2009) has shown the fraction of turbulent pressure to to- 
tal pressure on average increases up to ^ 30% at 2r5oo, 
roughly corresponding to rvir- They also found that the 
ratio at rsoo spans in a range from 6% to 15% for simu- 
lated relaxed clusters and from 9% to 24% for unrelaxed 
clusters, respectively. Dolag et al. (2005) has found that 
the turbulent energy accounts for ~ 5 - 30% of the ther- 
mal energy content. Vazza et al. (2009) has shown that 
the bulk energy changes from ^ 10% of the total energy 
at the center to ^ 30% at the virial radius and the tur- 
bulent energy ( 5%) is less than the kinetic one. The 
fractions of additional non-thermal pressure we derived 
above are slightly higher than those in simulated clusters. 
They should be considered as upper limits, taking into 
consideration the possibility of triaxial halo model, that 
decreases the amount of non-thermal pressure relative to 
total lensing one. 

We next computed entropy profiles from the lensing 
mass profile. Since the fraction of thermal pressure to the 
total one has not yet been revealed in reality, we cannot 
calculate the entropy of the thermal gas using lensing 
mass. We therefore compared the normalized gas entropy 
with lensing entropy (A'ions(f)), as shown in the right 
panel of figure 12. The normalization is applied at the 
third radial bin. The lensing entropy surprisingly agrees 
with the normalized gas entropy in the Offsetl direction 
out to the virial radius. The normalized entropies in the 
Offset234 and all regions also trace the lensing entropy 
up to '~ r^ooi while, in the outskirts, observed entropies 
deviate downward. Although we cannot constrain from 
figure 12 how much thermal energy is converted from the 
shock energy, when we remember that the slope of gas 
entropy in Offsetl also coincides with the accretion shock 
models (Tozzi & Norman 2001; Voit et al. 2002; Borgani 
et al. 2005, §4.4), our result indicates that the entropy 
generation by accretion shocks would occur, tracing the 
cluster mass profile. 

6. SUMMARY 

We observed Abell 1689 with Suzaku and detected the 
ICM emission out to the virial radius. Then, we derived 
radial profiles of temperature, electron density, and en- 
tropy, and conducted a joint X-ray and lensing analysis. 
Our conclusions in this work can be summarized as be- 
low. 
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The ICM emission out to the virial radius was de- 
tected with 4.0(7 significance. 

The temperature gradually decreases toward the 
virial radius from ^ 10 keV to ~ 2 keV. 

The temperature variations in azimuthal direction 
becomes significant in rsoo ^ r ^ rvh-: the tem- 
peratures of Offset 1 (- 5.4 keV) and Offset234 
(~ 1.7 keV) are about half and one-sixth of 10.4 
keV which is an emission-weighted temperature of 
2.'0-4.'0 and i'.O-6'.O regions, respectively. 

The temperature profile in the Offset 1 direction is 
in good agreement with the scaled temperature one 
by Pratt et al. (2007), while the temperatures in 
the other directions (Offset2 - Offset4) around the 
virial radius are lower than the scaled value by a 
factor of ~ 2. 

The electron density profile decreases down to 
0(10"^) cm~^ around the virial radius. 

The electron densities in the four offset directions 
are consistent with one another within errors in any 
circular annulus. 

The slope of the electron number density is — 2.29± 
0.18 in intermediate range of ?'25oo ^ r ^ f5oo- In 
the outer annulus of r^oo ^ r ^ Tvir, the density 
slope, — 1.24^0 52 where the lower limit is the sys- 
tematic error due to the contamination of X-ray 
emissions from point sources and inner region, be- 
comes shallower. This outskirts value is signifi- 
cantly shallower than the universal matter density 
profile p oc by Navarro, Frenk, & White (1996; 
1997). 

The entropies in all directions increase from the 
center to the radius of lO.'O. The entropy distribu- 
tion in the last radial bin (10.'0-18.'0) is anisotropic, 
similarly to the temperature profiles. 

The entropy profile of the Offsetl direction in- 
creases as r^--^, as predicted by models of accretion 
shock heating (Tozzi & Norman 2001; Ponman et 
al. 2003), while the entropies of Offset2-Offset4 in 
the outskirts are lower than those in their annuli of 
6.'0-10f0. The slopes of the entropy profiles in the 
Offsetl direction and the other directions match 
those predicted by lensing spherical masses with 
the hydrostatic equilibrium assumption, from the 
core to the virial radius and to ~ r^ooi respectively. 

The cumulative hydrostatic mass in the Offsetl di- 
rection increases with radius, while the ones in the 
other (Offset234) regions unphysically decrease out- 
side ^ 9.'0. The ICM in the cluster outskirts, ex- 
cept the Offsetl direction, is far from hydrostatic 
equilibrium. 

The anisotropic distributions of gas temperature 
and entropy in the outskirts are clearly corre- 
lated with the large-scale structure associated with 
Abell 1689. In the Offsetl (northeast) direction, 
the hot and high entropy region around the virial 



radius is connected to the overdensity region of 
galaxies in the SDSS map, whereas the low-density 
region (void) outside the virial radius is found in 
the other directions. It indicates that the ther- 
malization in the outskirts along the filamentary 
structure takes place faster than that in the void 
region. 

• The hydrostatic to spherical lensing mass ra- 
tio, MH.E.Z-^iens, Sit intermediate radii of 
3.'0 < r < 7.'0, is in a range of - 60 - 90%. 
Using the triaxial halo model, the X-ray to lensing 
mass ratio in the same annulus is consistent 
with unity within errors. To the contrary, the 
ratio within 2.'0, AfH.E.Z-Mions ~ 60%, is signifi- 
cantly biased low, and the discrepancy cannot be 
compromised by either the spherical or triaxial 
models 

• The gas fraction computed from combined gas den- 
sity and spherical lensing mass does not reach the 
cosmic mean baryon fraction within the virial ra- 
dius. Although our fraction is lower than the mean 
value of previous statistical studies based on lensing 
mass, the upper bounds of the gas fraction derived 
from the triaxial mass model (Oguri et al. 2005) 
come closer to, and are marginally consistent with, 
the statistical works. 

• The thermal gas pressure within ~ rsoo is at most 
~ 40 - 60% of the total pressure predicted by spher- 
ical lensing mass model. In the cluster outskirts 
outside rsoo; the thermal pressure in the Offsetl 
direction is comparable to total pressure, whereas 
it is 30 - 40% of the total pressure in the other 
directions. Taking into account the triaxial model, 
these values give the lower limits. With such val- 
ues the thermal pressure would be insufficient to 
balance the total gravity of the cluster, requiring 
additional pressure support (s). 
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APPENDIX 

SIMULATION OF SURFACE BRIGHTNESS PROFILE 

We simulated ICM emissions from Abell 1689 on XIS-FI (XISO + XIS3) with xissim, which spatially follow the 
doublc-bcta model determined from the 0.5-10 keV XMM-Newton MOSl image (§ 2.2). As a spectrum model, we 
used 6.2 keV apec model (metal abundance is 0.3 solar) with a flux of 2.5 x 10~^^ ergs cm^ s^^ (integrated within 30' 
in radius), which was determined from simultaneous fit to XIS-FI (averaged over XISO and XIS3) and XISl spectra of 
the Offsetl observation, extracted from whole XIS regions. Since Offsetl observation does not cover the central ^ 4' 
region, the spectrum model reproduces the emission from cluster outskirts. 

As shown in figure 13, the simulated profile is in good agreement with the observation. This indicates that the 
observed signals around the virial radius can be described by extrapolating the double-beta profile. From another 
simulation which used the same double-beta model but truncated at 10', we found that the stray signals which 
originate from < 10' regions contribute only 22% of signals in > 10' regions. 

In order to study an effect of temperature gradient expected in the ICM, we also simulated XIS-FI (XISO + XIS3) 
profiles of 9 keV and 2 keV thermal emission with apec model (metal abundance is 0.3 solar). When normalization of 
these two simulated profiles are scaled to match at the center, the 9 keV profile is higher than the 2 keV one by 10% 
outside 10'. Therefore, the temperature gradient does not affect the simulated profile significantly. 

SIMULATION OF POINT SOURCES 

Out of the 49 point sources in table 9, 35 sources have counterparts in the Incremental Second XMM-Newton 
serendipitous survey catalog (2XMMi, Watson 2009), on which fiux for 5 energy bands, 0.2-0.5 keV, 0.5-1.0 keV, 
1.0-2.0 keV, 2.0-4.5 keV,and 4.5-12.0 keV, are available from a HEASARC website^^. For sources which arc not on 
the 2XMMia catalog, we estimated their 0.2-2.0 keV flux by comparing their background-subtracted counts in the 
circular 1' regions in the PSPC total band (0.1-2.0 keV) image and that of source 7 (table 9) of which flux is known. 
From the estimated 0.2-2.0 keV flux, we further estimated flux for the 5 energy bands, assuming a power law spectrum 
with photon index of 1.4. 

Using the derived crude 5-band spectra of the 49 point sources, we simulated XIS-FI (XISO + XIS3) images for the 4 
observations using xissim. In order to improve statistics, we set exposures which are 10 times longer than the actual 
observations. As shown in flgure 14, the whole area of XIS is contaminated with signals from point sources. Figure 15 
shows the simulated 0.5-2 keV and 2-10 keV profiles of residual point-source signals in the four pointing directions. 

PARAMETRIC DENSITY AND TEMPERATURE PROFILES 

We parametrized the electron density and temperature profiles shown in figure 6, which is necessary to measure the 
hydrostatic equilibrium mass (§5.1) as well as to conduct a joint analysis of SuzakuX-ray satellite and Subaru/ Suprime- 
Cam+ HST/ ACS analysis (§5.2). 

The /3-model (e.g. Cavaliere & Fusco-Femiano 1976) or its generalized version (Vikhlinin et al. 2009) are often used 
to express electron density profiles in clusters. Since the low angular resolution of the Suzaku satellite does not require 
the complicated modeling (Vikhlinin et al. 2009), the analytic expression we use for deprojected electron density profile 
is given by 




(CI) 

where n^fi is the central electron number density, /? is the slope parameter, and r^^i and rc_2 are the core radii. In this 
modified expression of /? model, the number density of the outskirts r > lOfO are much better fitted with the outer 
radius rc,2, than a single /? model. Unfortunately, since the number of the data points is small, the core radius, rc^2, 
is not constrained well. We therefore fixed rc,2 = 3rc,i. Our results (§5) are not significantly changed by choosing a 
constant factor of rc,2/''c,i- The best-fit parameters for all, Offsetl and Offset234 regions are summarized in table 7. 
The observed projected temperature profile is described using a simple parametrization 
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with the temperature normalization Tq, the scale radius rt and the slope parameter a. We note that the assumed 
inner slope 1/2 is insignificant for our main purpose to describe the temperature profile in the outskirts. 

We fit the density and temperature profiles in the whole region of the cluster, Offsetl and Offset234 regions with 
the above two models, respectively. In the fitting of Offsetl region, we utilize the data of the other region instead of 
lack of the first two bins in the profiles. The best-fit parameters for all, Offsetl and Offset234 regions are summarized 
in table 7. 
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Fig. 1. — (left) A Mosaic image of Abcll 1689, created with 0.5 - 10.0 keV XIS-FI (XISO + XIS3) images smoothed with a Gaussian 
kernel of cr = 10". Background is not subtracted and vignetting is not corrected. CCD regions where *^Fe calibration sources are irradiated 
(Koyama et al. (2007)) are excluded, (right) the same as the left panel, but point-source regions arc specified as circles of 1' radius. See 
text for details of the point-source detection. 
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Fig. 2. — NXB-subtracted XIS spectra of N4636_GAL and Q1334, fitted with the X-ray background model described in § 3.2. Black and 
red are 1.5-7.0 keV XIS-FI (averaged over XISO and XIS3) and 1.5-7.0 keV XIS-BI spectra of N4636.GAL, respectively. Green and blue 
are 0.5-7.0 keV XIS-FI and 0.3-7.0 keV XIS-BI spectra of Q1334, respectively 





Fig. 3. — (left) Surface brightness profiles made of 0.5—2 keV XIS-FI (averaged over XISO and XIS3) image, using all the Offsetl through 
Offset4 observations. Point sources specified in figure 1 (right) are removed, but vignetting is not corrected. Raw (background inclusive) 
profile (black) is shown with XRB (green), NXB (magenta), and XRB+NXB (blue) profiles. A resulting background-subtracted profile 
(black — blue) is shown in red. The error bars in these profiles arc la. In the error bars in the red profile, Icr uncertainties of the XRB 
and NXB are added in quadrature to the corresponding statistical Icr errors, (right) the same as the left panel, but for 2-10 kcV profiles. 
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Fig. 4. — (left) A simulated 0.5-2 kcV profile of residual point source signals (gray), shown with background-subtracted 0.5-2 keV surface 
brightness profile of XIS-FI averaged over XISO and XIS3 (black). In the error bars in this black profile, Itr uncertainties of the XRB and 
NXB are added in quadrature to the corresponding statistical Icr errors, (right) the same as the left panel, but for 2-10 keV profiles. 




Fig. 5.— NXB-subtracted spectra of (a) 0'-2'. (b) 2'-4', (c) 4' -6', (d) 6'-10'. and (e) 10'-18' regions, which are extracted from Offsetl 
XIS-FI (black), Offsetl XIS-BI (red), Offset2 XIS-FI (green), Offset2 XIS-BI (blue). Offsets XIS-FI (cyan). Offsets XIS-BI (magenta), 
Offset4 XIS-FI (yellow), and Offset4 XIS-BI (orange). These spectra are fitted simultaneously with phabs X apec model, added to the XRB 
(GFE and CXB) model (§ 3.2). Although blank-sky spectra are also fitted simultaneously with the XRB model, they are not plotted here 
for clarity. See text for details of the spectral analysis (§ 4.3). 
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Fig. 6. — (a) Projected temperature profiles obtained by spectral analyses of Offsetl (red), Offset2 (green), Offsets (blue), and Offset4 
(cyan) observations. A temperature profile when all spectra of azimuthal regions at the same distance are summed is also shown in black. 
Dashed line is a temperature profile, kT = 12.3 — 7.7 {r/15'.6) keV, expected from the scaled temperature profile by Pratt et al. (2007). 
Dotted line shows the virial radius (15'6). Arrows are the positions of rsoo (8.'7) and r25oo (4'7). (b) the same as panel (a), but for 
deprojccted electron number density profiles. See text for the detailed method of derivation (§ 4.3). Dashed line shows a model profile by 
Peng et al. (2009) (Model 3 in table 9). (c) the same as panel (a), but for entropy profiles obtained by calculating K = kTjr?^^ from 
profiles in panel (a) and (b). Dashed line shows K oc r^'^. Normalization of dashed line is arbitrarily scaled. 




Fig. 7. — (a) Raw XIS-FI (black; averaged over XISO and XIS3) and XIS-BI (red) spectra of the Offsetl 10f0-18f0 region, and corresponding 
background spectra of XIS-FI (green) and XIS-BI (blue). Each background spectrum was created by adding NXB and XRB spectra. The 
XRB spectrum was simulated using fakeit command in XSPEC12, incorporating the GFE and CXB models detemined by the spectral 
fitting described in § 4.3. (b) The same as panel (a), but for the OfTset2 10.'0-18f0 region. 
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Fig. 8. — The hydrostatic equilibrium mass profile. The black, red and light blue colors represent the hydrostatic equilibrium mass for 
all, Offsotl and Offset234 regions, respectively. The solid and dashed lines indicate the best fit values and the 68% CL uncertainty errors 
obtained from Monte Carlo simulations, respectively. The green points are the model-independent mass based on the Chandra observation 
by Peng et al. (2009), which agrees with our results. The unphysical declines in the cumulative masses (all and OfTset234) indicate that 
the ICM of the outskirts is far from hydrostatic equilibrium we assumed. 
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Fig. 9. — The two-dimensional density distributions of galaxies around A1689, retrieved from Si? 55 photometric data, (left) : Distributions 
of red sequence galaxies of AbcU 1689. The box size is 40f0 X 40f0. The annulus is the outskirts region for spectral fitting, where inner 
and outer circles are lOfO and 18.'0, respectively. The central distribution is clearly elongated along the north-south direction. We could 
not find a significant difference between projected distributions in the Offsetl and the other directions, (right) : Distributions of galaxies 
whose photometric rodshifts span in the range of |z — Zc| < 0.016, where Zc is the cluster redshift. The boxes represent the FOVs of XIS 
pointings. The circle is the virial radius determined by lensing analysis. A Gaussian smoothing scale is represented by the hatched circle 
at bottom right in each plot. 
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Fig. 10. — Gas temperature map of Abell 1689 (shown in the black box), embedded in the 3 DSS gaXaxy map which is a magnification of 
figure 9 (right). Contours of the galaxy map are also overplotted for comparison. The temperatures are the best-fit values determined by 
the spectral analysis in § 4.3. See figure 6(a) or table 5 for their errors. 
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Fig. 11. — (left); Comparison of lensing mass, obtained from joint strong and weak lensing analysis, and hydrostatic mass. Blue, magenta 
and orange colors represent deprojected spherical lensing mass, NFW one, and triaxial halo model (Oguri et al. 2005), respectively. Black 
line represents the hydrostatic mass in all directions within the radius r ~ 7.'0, outside which the hydrostatic equilibrium is not adequate any 
more (figure 8). The solid and dashed lines arc the best fit values and the 68% CL uncertainty errors obtained from Monte Carlo simulations, 
respectively. Green points are the hydrostatic mass using Chandra data (Peng et al. 2009). The arrows denote the overdensity radii r2500 
and rsoo determined by deprojected spherical lensing mass, (right): Cumulative gas mass fraction. Blue, magenta and orange colors 
represent gas mass fractions using deprojected spherical lensing mass, NFW one, and triaxial halo model (Oguri et al. 2005), respectively. 
The solid and dashed lines are the same as left panel. Red points are the gas mass fraction based on AMiBA SZE measurement and Subaru 
lensing analysis for 4 clusters (Umetsu et al. 2009a), at three ovcrdensities A = 2500, 500 & 200. Light blue points are the gas mass fraction 
based on XMM-Newton X-ray and Subaru lensing analysis for 12 clusters of LoCuSS project (Zhang et al. 2010), at three ovcrdensities 
A = 2500, 1000 & 500. Green point denotes the gas fraction using Chandra hydrostatic mass for A1689 (Peng et al. 2009). In this figure, 
we use /i7i = Ho/71 km s~^ Mpc~^ = 1.0. 
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Fig. 12. — (left): Comparison of the gas thermal pressure with the total pressure expected from lensing masses with hydrostatic equilibrium 
assumption. Blue and magenta colors are the pressure calculated from deprojected spherical lensing mass and NFW one, respectively. Solid 
and dashed lines represent the best and ztlcr errors, respectively. Black, red and light blue points are observed thermal pressures (68% 
confidence errors) in the all, Offsctl and Offset234, respectively. Green line in the left panel is the pressure of Chandra data (Peng et al. 
2009). (right): The normalized entropy profiles. The colors are the same as those in the left panel. The slope of the entropy in Offsctl is 
in good agreement with that expected from lensing masses with hydrostatic equilibrium assumption. 




Fig. 13. — Simulated 6.2 keV apec (gray) profile is shown, which follows the double-beta profile determined from the 0.5-10 keV XMM- 
Newton MOSl image as described in § 2.2. For comparison, background-subtracted 0.5—10 keV surface brightness profile of XIS-FI averaged 
over XISO and XIS3 (black) is also shown. In this black profile, the residual contribution of point sources estimated in § 4.2 is also subtracted. 
In the error bars of the red profile, Icr uncertainties of the XRB and NXB are added in quadrature to the corresponding statistical la 
errors. 
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Fig. 15. — (top) Simulated 0.5—2 keV profiles of residual point source signals (gray) in OfTsetl through OfTset4 pointings from left to right, 
shown with background-subtracted 0.5-2 keV surface brightness profile of XIS-FI averaged over XISO and XIS3 (black). In the error bars 
in this black profile, la uncertainties of the XRB and NXB are added in quadrature to the corresponding statistical la errors, (bottom) 
the same as the top panels, but for 2-10 keV profiles. 
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TABLE 1 

Suzaku Observation log of Abell 1689. 



Start'' End=» 



Observation (ID) 



Oflfsetl (803024010) 
Ofrsct2 (803025010) 
OflFsctS (803026010) 
Ofrset4 (803027010) 



2008/July/23 04:42:40 
2008/July/24 08:23:04 
2008/July/25 10:09:59 
2008/July/26 11:24:00 



2008/July/24 08: 
2008/July/25 10: 
2008/July/26 11: 
2008/July/27 12: 



Exposure'' 



22 


23 


39.2 


09 


18 


38.2 


23 


19 


39.3 


13 


14 


37.9 



^ Time is shown in UT. 

^ Effective XIS exposure in units of ks, obtained after the event screenings. 



TABLE 2 

RASS AVERAGED COUNT RATES AROUND AbELL 1689 AND THE BLANK-SKY FIELDS. 



Band'' 


0.5°-1.0°'' 


Q1334<= 


N4636.GAL<= 


1/4 keV 
3/4 keV 
1.5 keV 


858.3 ±36.0 
157.6 ± 15.8 
139.1 ± 14.3 


816.1 ± 28.0 
133.3 ± 12.0 

132.2 ± 12.0 


1072.1 ±42.1 
238.3 ± 20.1 
159.2 ± 17.8 



"RASS 1/4 keV, 3/4 keV, and 1.5 keV bands correspond to 0.12-0.284 keV, 0.47-1.21 keV, and 0.76-2.04 keV, respectively 
'^RASS count rate in an annulus around Abell 1689, in units of 10~® counts s"'^ arcmin"^. 
•^RASS count rate averaged over the HLD field (1° radius) in units of 10^^ counts s^' arcmin^^. 



TABLE 3 

Results of model fittings to XIS spectra of the blank sky fields. 



r(cxB) 


5x(CXB)'' 


Sx(GFE)'' 


xVdof 


1.412 (fix) 


6.3+o;2 X 10"** 


3.5+i;2 X IQ-** 


162.4/127 



''The 2.0-10.0 keV CXB surface brightness in units of erg s ' cm ^ sr '. 
''The 0.3-1.0 keV GEE surface brightness in units of erg s^' cm^^ sr^'^. 



TABLE 4 

Significance of the signal outside 10'. 





Er action'' 


Significance'' 


Offsctl 0.5-2 keV 


0.21 


5.4o- 


Offset2 0.5-2 keV 


0.25 


2.9o- 


Offsets 0.5-2 keV 


0.45 


2.0fT 


Offsct4 0.5-2 kcV 


0.09 


2.6fT 


Offsctl 2-10 keV 


0.47 


2.1o- 


Offsct2 2-10 keV 


0.80 


O.Oo- 


Offsets 2-10 keV 


> 1 


O.Ocr 


Offset4 2-10 keV 


0.64 


0.2(T 



"Eraction of the residual point source signal defined as (simulated signal of point sources)/ (background-subtracted signal) outside 10'. 
''Significance of signal outside 10', after subtracting contributions from point sources and inter ICM emission. As to the latter contribution, 
constant value of 22% is used for all the pointings and energy ranges. 
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TABLE 5 

Fitting results of spectral analysis for the five annular regions. 



Region 

0'-2' 
2'-4' 
4'-6' 
6'-10' 
10'-18' 



8.42 
5.96 
5.37 



+2.76 
2.35 
.+2.89 
1.39 
+4.02 
-1.54 



0.77^ 



fO.OO 
' -0.77 

0.71+^" 



0.64 



-0.71 
1.07 
-0.64 



o cc + l.Ol 
0.00_Q 70 

8.39^111 

7 CI +3.40 
' ■'^^-2.00 

^■-'-^-o.so 

9 91+1.13 
^■^-'--0.83 



Z2" 



0.36^ 



f 0.15 
-'-0.15 

^2+°■^^ 

U.l^_0.i2 

o 7fi+0.00 

1.40+2-30 



0.49 



-1.17 
1.94 
-0.49 



11.391^7^ 
7.961^:68 

f 1.58 
-0.83 

1 3n+°-26 

l.JU_Q 24 



0.34"' 



^+0.09 
-0.09 
Q .^y+O.OO 



0.92"' 



4.49^ 



-0.77 
1+0.00 
-0.92 
-I 9-1+1.65 

o-5otJ:|| 



8.23 
10.53 
11.19 
6.27 
1.71 



,+0.46 
-0.42 
+2.60 
1.59 
+8.15 
■3.11 
+ 1.88 
1.21 
+2.56 
0.65 



0.26^ 



r. + 0.06 

-0.06 
49+0-26 
^•^^-0.25 

3.65l«:°° 



X^/D.O.F*^ 

752.0/605 
273.2/204 
319.0/210 
394.3/262 
665.7/458 



Region 



0'-2' 


='■■-''-'-0.39 


31+0-0-'' 
'^••3^^-0.05 


822.5/611 


2'-4' 


10.58^1-24 


63+0-30 


499.2/210 


4'-6' 


8.09ll:|g 


6^5+0-66 
U.OO_o.60 


429.7/219 


6'-10' 


c 01 +1-00 
'-'•■-'^-0.61 


U.9U_Q 4g 


423.0/271 


10'-18' 


9 Q9+0-85 
^•■3^-1.13 


^ =+0.65 
^•'J'^-0.31 


743.6/467 



^ fcTi , fcT2, fcTs, fcT4, and fcTaii represent for temperatures of Offsetl, Offset2, Offsets, Offset4 and the sum of all the azimuthal 
regions, respectively, in units of keV. Similarly, Zi. Z2 , Z^, Z^. and Z^\\ represent for metal abundances, in units of solar abundance 
of Anders & Grevesse (1989). 



Reduced x of the fit when parameters are separately determined for the four offset pointings. 
Reduced of the fit when parameters are tied between the four offset pointings. 



TABLE 6 

Deprojected electron number densities of the five annular regions. 
Region n^,!^ ne^2^ nc,3^ n.c,4^ 



'T-o,aII 



0'-2' - - 4.5310-06 X 10"' 4.7010-07 X 10-=* 4.62^0-0? x 10"' 

2'-4' - - 9.2310-24 X 10-4 1.0010-04 X 10-=* 9.6410-^2 x 10-^ 

4'-6' 3.0210;^;^ x IQ-* 2.77l0;2i x IQ-* 2.57l0::^4 x IQ-* 3.14l0;44 x IQ-* 2.88lO;33 x 10-* 

6'-10' 9.2ll};^^ x 10-5 8.44l2lg x lO"'"^ 9.86I2.I3 x lO^^ l-04lo;22 x 10"'' 9.4912.22 x IQ-^ 

10'-18' 5.1810-59 X 10-5 4.22li-^J x lO^'^ 4.2712-OJ x lO^^ 4.99li-2^ x 10-^ 4.66I1-33 x W'^ 

^ rici, nc,2, J^e,3, ^e,4, and ne,all represent for electron number densities of Offsetl, Offset2, Offsets, Offset4, and 
the sum of all the azimuthal regions, respectively, in units of cm^-^. 



TABLE 7 

Best-fit parameters for electron density and temperature profiles. 



Region 


n-o.o" 




rc,i" To^ 


e f 

rt 


All 

Offsetl 
Offset234 


(0.98 ±0.04) X 10-2 
(1.01 ±0.11) X 10-2 
(0.97 ±0.04) X 10-2 


1.16 ±0.02 
1.13 ±0.03 

1.17 ±0.02 


1.43 ±0.06 19.39 ±0.61 
1.36 ±0.13 15.85 ±1.30 

1.44 ±0.06 21.10 ±1.32 


1.03 ±0.07 3.88 ±0.29 
0.60 ±0.10 2.42 ±0.53 
1.32 ±0.19 4.72 ±0.69 


Ue^o, /3 and Vc^i are the central electron density (cm ^), the slope, the core radius (arcmin) for the 
parametrized density profile (eq. CI), respectively 

To, a, rt are the normalized temperature (keV), the slope, the scale radius (arcmin) for the parametrized 



temperature profile (eq. C2), respectively. 



TABLE 8 
Lensing Mass Parameters 



Model 




r - b 

wvir 


NFW^ 
Spherical"^ 
NFW-equivalent triaxial model" 


1 47+0-59 

4-^'-0-33 

1.0U_Q 3g 

1 26"'"0-28 
4-^0-0.49 


12.7 ± 2.9 

16.9l?2'9 



^Virial mass. 

Concentration parameter of the NFW profile. 
■^NFW model parameters taken from Umetsu & Broadhurst (2008). The quoted errors include both statistical and systematic uncertainties. 
"^Deprojected spherical mass. 

"^NFW-equivalent triaxial model parameters taken from Oguri et al. (2005). 
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TABLE 9 

Point sources identified in XIS FOVs of OffsetI through Offset4 observations. 



4J.a 

rr 








Flux*^ 


1 * 


lyo.uo^z / 


1 
- 1 


1 91 Q1 p. A 
. IZIOID^L 


Q OOzlOf^r. 1 A 
O.Uy4UOO- 14 


A 


1 OQ nccc; 1 /I 

lyo.UOol4 


- 1, 


t^K0997fi 

.DDyzz i 


OAA GA 1 A 
Z . 44y4De- 1 4 


o* 
O 


lyo.uzoyo 


1 

- 1, 


.ozooy 1 Z 


9 t^7QQQa 1 A 

z.o / ooot;- 14 


A 


ly i .ozo'y:o 


1 

- 1, 


941 R7C\^ 
.Z^IO / uo 


1 1 OOOQo 1 

1 . ivjuuyti- lo 


O 


1 Q7 7^077 


1 

- 1, 


.^oZoO'y: ( 


7 70Rfi<-> 1 A 
1 . 1 UDDC-14 


a 
D 


ly i . i y^oo 


1 

-1, 


c;n9'^09Q 

.ouzovjzy 


Q CI 70^0 1 A 
O.Ol 1 ODti- 14 


7 


ly 1 .yoD^o 


1 

- 1 , 


990*^008 

.zzuoyuo 


A 71 SAzlo 1 A 
4. 1 lo440-14 


Q 

o 


ly 1 . i Doiz 


1 

- 1 , 
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44 


197.81756 


-1, 


.1859025 


2.09468e-14 
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^Serial number for point sources. Sources with * on the number do not have counterparts in the 2XMMi catalog. 
^Positions of point sources in the Equatorial ,12000 coordinates {degrees). 
'^The 0.2-12.0 keV fiux in units of erg s"^ cm~^. 



